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The separation of the Weyl nodes in a broken time-reversal symmetric Weyl semimetal leads to 
helical quasi-particle excitations at the Weyl nodes, which, when coupled with overall spin conserva¬ 
tion allows only inter-nodal transport at the junction of the Weyl semimetal with a superconductor. 
This leads to an unusual periodic oscillation in the Josephson current as a function of koL, where 
L is the length of the Weyl semimetal and 2ko is the inter-nodal distance. This oscillation is robust 
and should be experimentally measurable, providing a direct path to confirming the existence of 
chiral nodes in the Weyl semimetal. 

PACS numbers: 74.45.+c, 74.50.+r, 73.21.-b 


Introduction .—Weyl semimetals (WSM), which have 
received much interest recently due to their non-trivial 
transport characteristics, are 3D topological systems 
where conduction and valence bands touch at two or more 
‘Weyl’ point^iM^. According to a no-go theorenP, gapless 
Weyl nodes in a WSM appear as pairs in momentum 
space with each of the nodes having a definite ‘chiral¬ 
ity’, a quantum number that depends on the Berry flux 
enclosed by a closed surface around the node. Gauss 
law prevents the annihilation of the nodes unless two of 
them with opposite chirality are brought together, which 
provides the ‘topological’ protection of the Weyl node^. 
A WSM phase requires broken time-reversal and/or in¬ 
version symmetry and a growing number of systems has 
been put forward which realize the WSM phas^^®^. 

The separation of the chiral nodes, allows charge 
pumping between the nodes in the presence of parallel 
electric and magnetic helds, as a consequence of the chiral 
anomaljffil, and this has led to detailed studies of trans¬ 
port in Weyl semi-metals in several recent paperJi^l^. 

In this paper we study the current in a simple Joseph¬ 
son junction setup, depicted in Fig. Ha). The heli¬ 
cal quasi-particle excitations at the Weyl nodes, due 
to the overall spin conserving processes at a WSM- 
superconductor (SC) junction, allow only inter-nodal 
transport!^. Further, we show, unlike in a normal metal- 
SC interface, the inter-nodal ‘normal’ (electron to elec¬ 
tron) reflection process in a WSM-SC interface is not 
suppressed even for energies close to Fermi-energy, due 
to the broken time-reversal symmetry separating the 
Weyl nodes. The Josephson current, flowing through the 
bound levels formed by multiple inter-nodal ‘normal’ and 
Andreev (electron to hole) processes in a SC-WSM-SC 
system, consequently, acquires a specific periodicity as a 
function of the length of the WSM which depends only on 
the separation of the Weyl nodes in the momentum space 
(see Fig. Hb)). We argue that both of these features are 
robust because they are not only bulk effects, but they 
are also protected by the robustness of the Weyl nodes. 
We also discuss the feasibility of experimental observa¬ 
tions of this transition in our system, which can confirm 
the presence of chiral nodes in WSM. 



FIG. 1. (Color online) (a) The setup for Josephson current 
with two superconductors (SC) characterized by phases (pi, 
and (pn sandwiching a WSM of length L between them where 
a time-reversal broken perturbation separates the Weyl nodes 
in momentum space by 2ko in k^. (b) The Josephson current 
(at normal incidence) is periodic in L with a period of nfko. 
We show the zero-temperature Josephson current, Eq. §, 
as a function of the superconducting phase difference (p, for 
various values of L in solid (dashed) lines for 9 between 0 and 
TT (tt and 2tt), where 9 = 2koL mod(27r). The parameters used 
are koL = 31.4,32.0,32.5,33,33.5,34.0,34.454, I?kll2mw = 
lOpiv = lO^A = p-sl2., p <C feo and ms = mw- 


This oscillation in the Josephson current and the re¬ 
sulting changes of sign of the critical current at arbitrary 
values of (p (or the O-tt transition) is an inherent prop¬ 
erty of the SC-ferromagnet-SC junctiorP^^^^and has also 
been experimentally observecpSl. Since our model also 
explicitly violates time-reversal invariance, our results 
show quite a strong similarity with the Josephson current 
in similar systems^ as well in semiconductor nanowires 
with Zeeman couplinj^. 

Model and geometry .—We consider the geometry as 
shown in Fig. [^a) with the superconductors at z < 0 
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and z > L and the Weyl semimetal (WSM) in the region 
0 < z < L. We model the WSM starting from the stan- 
d ard H amiltonian describing a 3D TI in the Bi 2 Se 3 fam- 
il^ l4i | 42 | regularized on a simple cubic lattice and adding 
a time-reversal breaking perturbation bz to access the 
WSM phas^i^- 


Ho =ekTx - Xz sin kzTy 

- XTzieJxSinky - aySinkx) + bzaz- ( 1 ) 

Here ek = e — 2t cos ki is the kinetic energy, r {a) 
represent the orbital (spin) degrees of freedom and A, 
A^ are the strengths of the spin-orbit coupling. In the 
limit Az <C M <C fez, (where M is defined as e -ePH), 
this simplifies to a two-band model for a WSM, where, 
in the absence of t he sp in-orbit coupling. A, the bands 
have opposite spinP^HIl xhe chiral fermion excitations 
around the two Weyl nodes (which we choose to be at 
(0, 0,±fco), where tk^ = bz — M) are described by the 
Hamiltonian 


HwSM — ^k^z X{kx^x ky^y) ( 2 ) 


with Ck = {b?/2mw){kl + ky + k^ — k^) being the ki¬ 
netic energy, /iw being the chemical potential measured 
from the Weyl node and my/ being the effective mass. 
Since the Weyl nature of the fermions is only evident at 
momenta which are small with respect to the symmetry 
breaking scale fcoj we choose the Fermi energy ^yy to be 
small enough, so that the Fermi surfaces around the two 
Weyl nodes are disconnected. The surface states for this 
model appear on the surfaces perpendicular to the x-axis 
and y-axis. In this paper, we do not attempt to access 
transport through the surface states. Instead, we con¬ 
sider transport through the bulk of the WSM and hence 
restrict ourselves to transport parallel to the z-axis. 

The superconductors can be described in terms of the 
Boguliobov-de Gennes (BdG) Hamiltonian as 


Hic = 


^kl2x2 e^‘*'^Aiay\ 
-e-^'f’^Aiay -^khx2 j ’ 


( 3 ) 


where A is the pairing potential in the superconductor 
and £,k = (^^(^^ + ky-\- kl)/2ms — ns)- ins is the effec¬ 
tive mass of the electron in the superconductor and y,s is 
the chemical potential. is the superconducting phase 
of the superconductor. For the left and right super¬ 
conductors, j = L, R. The parameter depends on the 
details of the superconducting material. In the numerical 
results shown, we consider /xg ^ A, which is the realistic 
limit. Also, for simplicity, we consider mg « myy- 

WSM-SC junction .—The solutions of Eq. Q in the 
Nambu-Gor’kov space are now 4 component spinors. For 
incident energy E, the right-moving solutions with the 
wavefunctions proportional to for electrons and 

^-vik^z holes can be written in the basis of the two 
bands iz = ±, with 


K{h) = \/fco + vi2.m.ui/b?)\/{pw + (-)A)2 - (Ap)2 




FIG. 2. (Color online) (a) Both reflection (R) and Andreev 
reflection (AR) in WSM occur from one chiral node to an¬ 
other. The chiralities of the nodes are denoted as -fve and 
—ve, whereas the solid and the dashed lines show dispersions 
of Eq. with positive and negative velocities (= dE/dk). 
The bands have opposite spins, which accounts for the change 
of chirality for both ‘normal’ and Andreev reflection, (b) The 
probability that an electron will be reflected as an electron 
(Ix^P) or 8-s a hole (|?7^P) at a WSM-SC interface, dis¬ 
cussed in Eq. Note that the probability of reflection as 
an electron is finite. The parameters used are the same as in 
Fig.^b). 


and p = yk'j. + . The left-moving solutions can be 

written similarly with —>■ —For the case of a 

WSM-SG junction, the WSM and and the superconduct¬ 
ing wavefunctions on the two sides of the junction can be 
matched at the junction by requiring the continuity of 
the wavefunction and its first derivativ^^. This leads to 
the net reflection matrix W from the WSM-SC junction, 
which connects the left and right-moving solutions. 


= 




( 4 ) 


where the ‘normal’ reflection matrices and the 

Andreev reflection matrices ?'g^(r^g) denote, respectively, 
electron to electron (hole to hole) and hole to electron 
(electron to hole) processes at the interface with the jth 
superconductor. 
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FIG. 3. The variation of the bound levels (solutions of Eq. §) 
near the chemical potential with the length L of the WSM for 
various values of 9, where 9 = 2koL mod(27r). The parameters 
used are the same as in Fig. Sb). 


For the purpose of physical interpretation, let us take 
the case of near-normal incidence (fco 3> p) of an electron, 
where the reflection matrices reduce to the form: 


r 


•7 = 

ee 





(5) 


In this simplified form it is immediately clear that both 
the reflection and the Andreev reflection change the chi¬ 
rality (see also Fig.j^a)) and can only take place from one 
node to another because of the chiral nature of the nodes. 
We plot the probabilities of normal and the Andreev re¬ 
flection in Fig. I^b). We note that even at energies close 
to the Fermi energy, normal reflection is not suppressed. 
The existence of the new momentum scale fco 7^ in¬ 
troduced by breaking the time-reversal symmetry, allows 
the incident electron momentum to be different from the 
Fermi momentum of the supercondu ctor. This leads to 
the non vanishing of normal reflectiorP2l211 

In contrast, note that for a topological insulator in 
3D, the bulk is gapped and the non-trivial transport in 
junctions with superconductors is purely due to the sur¬ 
face states, where, the surface states consist of a Dirac 
metal with an odd number of nodes whose fermions have 
their spins aligned with the direction of motion (spin- 
momentum locking). This leads to completely different 
ph ysics f or a topological insulator-superconductor junc- 
tioiPUm Two dimensional graphene, on the other hand, 
is metallic and the transport is through the bulk. How¬ 
ever, in graphene, although there are two Dirac nodes 






FIG. 4. (Color online) The Josephson current as a fnnction of 
both L and fco is shown at the value of 0 « 7r/2. The initial 
valne at the origin is {ko,L) = {q,l), ql ~ IOtt. The contonrs 
of constant current follow a set of (approximate) hyperbolas 
for constant 9 = 2koL mod(27r), a few of which are shown 
in the right margin (with the minimum and the maximum 
current occuring near 9 = n and 9 = 0 respectively). Other 
parameters used are the same as in Fig. Sb)- 


(valleys), each of the nodes has fermions of both chiral¬ 
ities and the resulting process at the a superco nduct ing 
interface is purely intra-nodal Andreev reflect ioip2llll_ 

Bound levels in the SC-WSM-SC geometry .—Multiple 
reflections at the WSM-SC boundaries lead to bound 
electronic levels in the SC-WSM-SC geometry. But as 
discussed above, normal reflection amplitudes are not 
small at the WSM-SC interface and, in general, there 
is no simple way of summing up the amplitudes between 
the two superconductors to obtain the resonance con¬ 
dition when both Andreev and normal reflection ampli¬ 
tudes are non-zero. For the case of near normal incidence, 
however, the problem simplifies and the bound levels Eh 
can be found by solving 


det [/ 4 x 4 - ■R^MTZ’^M] = 0, (6) 


where M is the matrix which accounts for the phase the 
electron/hole acquires while moving from one junction to 
another. We note that for ko 3> p, Eq. ([^ can still be 
used for approximate solutions. Writing 


n^MTZ^M 


%e Teh 

The Thh 


(7) 


in the limit of near normal incidence with /cg I2mw much 
larger than incident energy E and ps much larger than 
pairing potential A, the T matrices have the simplified 
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form (with ms = my/)- 


T = 

• ee 


a+ 0 
0 a~ 


The = 


0 /3+ 
/ 3 - 0 


with « g± 2 ifeoi , 

/3± « ±e=F2*'=“-^2i(l + e"*'^)A, 


( 8 ) 


where ^ = y/2rnw~{^/4> = 4 >r — 4>Lt Also Thh = 
Tee{E —)■ —E), Teh = TheiE ~T)- This immediately 
shows the periodicities of the T matrices, T{(j)) = T{(j) ^ 
(j) + 2 tt) and T{2koL) = T(2koL —>■ 2fcoT + 27r), which 
implies that the bound levels Et, the solutions of Eq. 
also inherit the same periodicities in 4> and 2koL. This 
additional periodicity of the levels with period (7r/fco) 
in length appears as a consequence of the inter-nodal 
normal and the Andreev reflections. The periodicities of 
Eh in the difference of the superconducting phases (j) and 
in 2fcoT, in the limit of fco 3> p is shown in Fig. This 
is our central result. 

Periodic oscillations in the Josephson current. —The 
Josephson current for the system with the total Hamil¬ 
tonian El is written as Jjos = x (w)’ 
erage is taken over the states of the system. For the 
non-interacting system, where the length L is much 
smaller than the coherence length in superconductors, 
the Josephson current flows through the bound levels 
(neglecting the continuum contribution) and can be esti¬ 
mated as^ 

o „ n jT 

J{ew) =-^f{Eb - Ew), ( 9 ) 

where / is the Fermi-distribution function. Apart from 
the 27r periodicity of the Josephson current in (p, as the 
bound levels Eh are periodic in L with the periodicity of 
7r/fco, the Josephson current also inherits the same peri¬ 
odicity. This periodicity is shown explicitly in Fig. m 
for the case when p = 0 . 

The periodic dependence in L can also be written as 
an approximate periodicity in fcg with a period of tt/L. 
For large values of L, the rapid oscillations of Eh with a 
small variation of ko outweighs any other dependence on 
fco and the periodicity is almost exact. The Josephson 
current as a function of both /cq and L is also shown in 
Fig.a where the locus of constant current approximately 
follows 9 = 2koL mod(27r). This is another of our main 
results. 

Lattice simulation. —To analyse the case for non¬ 
normal incident angle p 7 ^ 0, we compute the Joseph¬ 
son current from the lattice version of Eq. Q through 
its Green’s functioij J^Fo i Green’s function of the 

WSM, g{oj) = [{uj + iS)I — E[o]~^, is coupled with two 
superconductors i = L,R through the on-site self en- 
ergj^ 


E,(w) = (10) 



FIG. 5. (Color online) Result of lattice simulation: (a) 
Josephson current as a function of the superconducting phase 
difference p for various values of 0 = 2koL mod(27r). (b) 

Josephson current as a function of the length (kpL/ n) at 
(f> = tt/2. The parameters used ard^t = 0.25t (see Eq. ( |10[ |), 
e = 6t, A = Az = t, A = O.Olt, fi — 0.02t and L — QO sites. 
For these parameters, the position of Weyl nodes, ko — bzjXz 
(the lattice constant is taken as the unit of length). 


where C acts on the particle-hole degree of freedom 
in the Nambu basis and is defined only on the 
sites in contact with the superconductor, t charac¬ 
terises the tunnelling between the superconductor and 
the WSIVP^. Then writing the full Green’s function as 
G{oj) = (p“^(w) — Er {u}) — T,r{uj))~^ we compute the 
Josephson curreniP^Ml Wg gncl that the oscillations of 
current remain intact and this result has been summa¬ 
rized in Fig. The bound levels Eh can also be found 
approximately for p ^ ko using Eq. , and the corre¬ 
sponding Josephson current also shows that the oscilla¬ 
tion with kgL remains intaclP^. 

Feasibility of experimental reali zation .—In the pre¬ 
dicted WSM material TaA^^U^, chiral node pairs 
(formed by breaking inversion symmetry) are separated 
in momentum space by a distance ~ 0.02A“^. Assuming 
standard electron mass, the relevant energy scale is about 
a milli electronvolt, which only becomes larger if the ef¬ 
fective mass is smaller. Gombining this with the fact that 
large momentum scattering (from — fcg to kg) is needed to 
break the topological protection of the chiral nodes, heli¬ 
cal excitations in WSM are expected to be robust against 
disorder in a relatively clean sample. The periodicity of 
the Josephson current as well as the bound levels that we 
have discussed are, in principle, observable in tunneling 
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experiments. The periodic variation of the bound-levels 
can also be probed in Andreev spectroscopy. For a typi¬ 
cal sample, the length scales for such periodic variations 
would be of the order of few tens of nanometers. The 
separation of the Weyl nodes can also be tuned by ad¬ 
justing the magnetic dopin^sHHH ^ observe periodicities 
with the separation of the Weyl nodes. 

The effect of having many Weyl nodes complicates the 
theoretical modeling and presents a weakness in our pro¬ 
posal. But, as long as the transport takes place along a 
pair of Weyl nodes, a similar periodicity in the Josephson 
current is expected. 

Summary and conclusion .—To summarize, we have 
shown explicitly, employing a simple model of the WSM, 
the occurrence of inter-nodal reflection processes at an 
WSM-SC interface due to spin conservation. This gives 
rise to an unusual periodicity in the bound state spectra 


and consequently in the Josephson current that depends 
only on the separation of the two Weyl nodes and the size 
of the sample. This provides a direct path for possible 
observations of the manifestation of inter-nodal Andreev 
reflection in Weyl semimetals. 

In closing, we sketch some problems for future studies. 
Apart from transport signatures of the chiral anomaly in 
WSM, the appearance of surface states, and the conse¬ 
quent Fermi arc dispersion is another remarkable feature 
of the time-reversal broken WSM. Their transport char¬ 
acteristics in the Josephson current would be interest¬ 
ing to study. Finally, quantitative investigations of the 
effects of disorder and interactions on transport in the 
WSM are also left for future studies. 

A. K. was supported in part by the NSF through Grant 
No. DMR-1350663. 
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APPENDIX 

A. Model hamiltonian for the WSM 

In this section we derive Eq. (2) near the Weyl nodes 
in Eq. (1) of the main text. Starting with parent Hamil¬ 
tonian Eq. (1): 

H{k) =ekTx - X^sink^Ty 

- Atj (tJa; sin fcy - CTj, sin fca;)-b (11) 

we define a mass term M = e — 6t, so that, around the E 
point 

Ck = M + 2t ^ ^ cos Ki M + t ^ kf. (12) 

Eor X = 0 = the hamiltonian (about E point) is (M + 
t\k\'^)Tx — XzkzTy, giving two doubly degenerate bands 
with the dispersion 

« ±^M‘^ + kl{2Mt + XI) + 2p2(Mt), 

where = k"^ + ky. Assuming that the mass is large 
enough, we can write 

± a/M 2 -b kl{2Mt + A2) -b 2p2(Mt) « ±{M + akl+tp'^), 

where a = t + X^/2M. This quadratic dispersion with a 
non-zero gap (2M) at the E point is the starting point in 
Ref. [33l Motivated by this, we define a unitary matrix U 
such that 

[{M + t\k\‘^)Tx - XzkzTy] U = {M + akl -b tp^)Tz. 
In this new basis the full hamiltonian WH{k)U is 


( M + bz + ak^ + tp^ 

I 

X{ky - ikx) 


0 

M — bz + akl + tp^ 

X{ky ~\~ ikx) 

0 


0 

X{ky ^kx) 

—M + bz — akl — 

0 


tp"^ 


X{ky -b ikx) 

0 

0 

—M — bz — ctk^ — 



(13) 


Clearly, bz lifts the spin degeneracy and shifts the bands 
up or down in energy while A mixes the uppermost and 
lowermost bands and the two bands in between. If the 
chemical potential is not large, the lowest energy exci¬ 
tations are only in the two bands close to zero energy. 
Limiting ourselves to the lowest energy excitations, we 
get the effective 2 -band model 


( t{a'kl + p^ - kl)w X{ky-ikx) \ 

X{ky+ikx) -t{a' kl + p^ - kl) ) ' 


where tk^ = bz — M (assuming bz > M) and a' = a/t = 
1 + Xl/2Mt « 1. This is the hamiltonian that we use in 
the main text. 


B. Solving the WSM-SC interface 

In this section, we provide the details for the derivation 
of the reflection matrix in a WSM-SC system. Eollowing 
Ref. 35 (Uchida et ah), the wavefunction of energy Ei in 
the WSM is given by the following solutions of Eq. (1) 
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of the main text in the Nambu-Gor’kov space (with the 
Hamiltonian in the hole space written as — 

i’wsuiz) = E I 
a=± I 

+'H'^ |, (15) 

where a is the band index, a{b) denotes the electron 
(hole) amplitude and L{R) denotes the left (right) mov¬ 
ing solution. are normalized eigenvectors, which 

are non-zero in electron (hole) sector of the Hamilto¬ 
nian. In each sector £{T-L)^ oc (/e(?i), (—)A+(_))^, and 
EiUy oc ((-)A_(+),/e(h))^, with /e(h) = ^lw + {-)Ei + 
V^(/rvv+7-pA)2^^0Ap)2, A± = \{k^ + iky). 

In the superconductor, the solutions of Eq. (2) of the 
main text are: 


C. Non-vanishing normal reflection 


For the simplest case, if we take the only incident 
amplitude to be nonzero as aj = 1 on the left side of 
Eq. (161, we can show that the amplitude of normal re¬ 
flection is 



with Md; = vc^ ^1 - / ^-1 - • (17) 

In a normal metal (with normal incidence), qg Ki 

K, kp when the incident energy is small, so that the 
amplitude of the normal reflection a~^ vanishes. However, 
in the WSM, q^ Pi q^i k kp, but kf Ki ko ^ kp. Hence, 
for the WSM, the normal reflection is not suppressed in 
general at small incident energies. 


fpsc{z) 


uc^ 

UCl 

-VCi 

vc^ 


( 

—vd^ 

ud^ 

V udi 




where, with £1 = ~ Ef, 

u{v) = + {-)m) /2E, 


D. Bound state spectrum 

The bound state spectrum for the SC-WSM-SC geome¬ 
try can be found by using the reflection matrix at a single 
WSM-SC interface defined in Eq (S6). At near normal 
incidence, the bound states are given by the zeroes of the 
determinant 


and qe and —qu are, respectively, the outgoing electron 
and hole momenta in the superconductor, defined as 
(with Fermi momentum kp) 

qe{h) = sjkp-p"^ + i-)2msi£l/k? . 

The boundary conditions at z = 0 are given by the 
continuity of the wavefunction and its derivative: 


det [/ 4 x 4 - n^MTZ^M] = 0, (18) 

The TZ^ and TZ^ matrices are the reflection matrices at 
SC-WSM and ITS'M-5'G interfaces respectively and the 
J\4 matrix is the phase picked up by the electrons and 
holes in the WSM region (of length L). At near normal 
incidence, we can compute the matrices analytically and 
write 


V'wsm(O) = tpsciO) 


ms 


Uz 0 
0 (Jz 


^zV'WSm(^) | 2=0 = 


with Gz being the Pauli matrix, 
gets the reflection matrices. 


( at 


\ 


\K I 


Tee Teh 
^he ^hh 


mwdztpsciz) \z=o, 

By solving them one 


/ 4 \ 

\^R / 


(16) 


TZ^MTZ^M 


%e Teh 

The Thh 


where {(j) 
Tee = 
The = 


ti? - 4‘l), 

at + 0 \ 

0 a)" -|- a^ J ’ 

0 ;3+(e-*^/3++/3+) 

Po +P 2 ) 0 


± _ ‘ie^dkt+k^)L^-i,Pi^±j^T(^q^ y qh)Vv^ 

[{ke + qe){k^ + qh)yP - (fc/t - qe){kt - qh)v‘^]^ 

^ [{kf - qe){kl E qh)u^ + {-k^ + qe){kf + qh)v^Y 

[(fct -f qe){k^ E qh)u^ - {k^^ - qe){ke - qh)v^] ^ 















/3^ = T‘2kf{qe + qh)uv, 

P± ^ -^i(kT + k^)L [(fc? + ge)(fc^ - qhW - (k^ + ge)(fcf - qhW] 
[{kf + qe){k^ + qhW “ {k^ “ 9e)(fc? - qhWY 
P± ^ ^^2rkfL [(^? - + <lh)v? " " ge)(fcg' + qh)v^] 

[(fcj + 9e)(/c^ + qhW - {kh - 9e)(fcf - qhW] ^ 



FIG. 6. (Color online) The zero-temperature Josephson cur¬ 
rent at a fixed superconducting phase difference ^ = -k 12 with 
the length of the WSM, L, keeping kg the same for transverse 
momenta fca;, fcj, = 0 in the blue (dashed) curve, while having 
the transverse momenta = 0.1, fcy = 0 in the red (solid) 
curve. Other parameters are the same as in Fig. 1(b) of the 
main text. 


and Thh — 'T'ee{k!j ~E), 'Teh — —E). This 

immediately shows the periodicities of the E matrices, 
E{4>) = T{4> —>■ (j) ~k 27r) and E(2koL) « T(2koL —>■ 
2ki)L + 27r) considering E <C k^/mw, he, k^^^^ « fcg- 
The expressions can be simplified in the limit /ig A 
and fcg/2mry E (the incident energy) and are given in 
the main text. 


E. Variation of Bound levels for p 7 ^ 0 

We show a plot of the zero-temperature Josephson cur¬ 
rent with varying length of the WSM in Fig. for trans¬ 
verse momentum p ^ 0 (but still one order of magnitude 
smaller than fcp). For the approximate determination of 
the bound levels Eb, we assume that Eq. ( 6 ) of the main 
text still remains valid in this range of parameters. The 


figure shows that the periodicity with k^L is intact, which 
implies the robustness of the periodicity in non-normal 
but small-angle incidence of electron. 


F. Green’s function formalism of Josephson current 


Following Ref.[^ if W represents the number operator 
of the ith site in the metallic system (WSM), then —eW 
can be written as the sum of the current flowing from the 
{i — l)th site to the ith site and from the fth site to the 
{i + l)th site, i.e, W = (—l/e)(+ Ji+i^i). Each 
of these terms represents the Josephson current flowing 
in the system and is independent of the site i for a large 
enough system. To elaborate, let us write the lattice 
Hamiltonian for the WSM in the following form (in the 
basis of spin operator 172 with eigenvalues a and parity 
operator Tz with eigenvalues r): 


H = H^ + H^= y at /i y a 

/ ^ ICTT t,crcr T,T 


,r,r' 


V a} a / 

(ij) ,a,a',T,T' 


Then first term does not contribute to the current. In 
our case the second term of the Hamiltonian is (assuming 
translation invariance in the directions perpendicular to 
that of the flow of Josephson current) 

H'^ = al{-twSMTx +iXzTy)(Toaj + h.c. . ( 20 ) 

And the current from the ith to the {i + l)th site as 


Ji{t) = y (^WSM — T Az)x 

£T,r 

where f = —r. The above averages can be written in 
terms of the Green’s function as 


Q+- 
^ij,a<7' ^TT 





t 

i(7T 

t 

i(TT 




(t)) (n<2(7'7-(^ 

(t)) {aiurit 


{t))\ 
(t)) ) ■ 


In the absence of any applied voltages the correlation 
function G+“(a;) is given by G+“(a;) = /(w)[Ga(w) — 


Gr(w)]. /(w) is the Fermi function and Ga/r are the 
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advanced and retarded Green’s functions. 



